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1 Introduction 

Fragmentation is natural phenomenon that can be observed at a great variety of scales. To 
give just a few examples, we may think of stellar fragments and meteoroids in Astrophysics, 
fractures and earthquakes in Geophysics, crushing in the mining industry, breaking of crystals in 
Crystallography, degradation of large polymer chains in Chemistry, fission of atoms in Nuclear 
Physics, fragmentation of a hard drive or files in Computer Science, ... In this text, we will be 
interested in situations where this phenomenon occurs randomly and repeatedly as time passes. 
Typically, we may imagine the evolution of blocks of mineral in a crusher. 

Over many years, the importance of the activity on fragmentation in Physics (see e.g. [21I2T] 
and the references therein) has not been reflected in Probability Theory where the interest 
has been more irregular. The first significant works concerning its probabilistic aspects (strong 
laws) are due to Kolmogorov j3H] and his student Filippov . The systematic study of general 
fragmentation processes is quite recent, although some well-established areas of Probability 
Theory, such as branching random walks and multiplicative cascades, are clearly relevant to 
investigate a large family of fragmentation processes. 
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This text is intended as a survey on a stochastic model for fragmentation processes which 
is characterized by a few parameters. We shall present various point of views on this model, 
from which different properties depending on the values of the parameters are derived. In order 
to deal with models that can be studied mathematically, we are led to make hypotheses that 
may look at first sight somewhat stringent, but which are often fulfilled in applications. First, 
we suppose that the system has a memoryless evolution, i.e. its future only depends on its 
present state and not on its past. In other words, the system enjoys the Markov property. In 
particular, this excludes the possibility that an object might be more fragile (i.e. more likely to 
split) due to former shocks. Second, we assume that each fragment can be characterized by a 
real number that can be thought of as its size. This impedes to consider geometrical properties 
like the shape of a fragment. Third, we shall suppose that the branching property is fulfilled, 
in the sense that fragments split independently one of the other, or in other words, that the 
evolution of a given fragment does not depend on its environment. Finally, we shall assume 
that the process enjoys self-similarity, that is that the law of a fragmentation process started 
from a unique fragment of size r > can be reduced by proper rescaling in space and time, to 
that when r = 1. 

In order to give a formal definition, we introduce the state space of decreasing numerical 
sequences which tend to 0, 

S := \ s = («!, s 2 , ■ . .) : Si > s 2 > . . . > and lim s n = > . 

A generic configuration s should be thought of as the ranked sequence of the sizes of some object 
that has been split. The space S is endowed with the uniform distance, so S is a complete 
metric separable space. 

Let X — (X(t),t > 0) denote a random process such that 

X(t) = (X 1 (t),X 2 (t), ...)eS Vt > 0, a.s. 

We suppose that X is continuous in probability and Markovian. For every r > 0, we write P r 
for its law started from the configuration (r, 0, . . .), and we further assume that for every t > 0, 
Pi(Xi(t) < 1) = 1 and ¥ (Xi(t) = 0) = 1. In other words, the sizes of the pieces resulting 
at time t from an object with initial size 1 cannot exceed 1, and an object with size cannot 
produce pieces with positive size. 

Definition 1 We call X a self-similar fragmentation if : 

(i) There exists a G R, called the index of self-similarity, such that for every r > 0, the 
distribution under Pi of the rescaled process (rX(r a t),t > 0) is P r . 

(ii) For every s = (s 1 ,s 2) ---) £ S, if (X^ l \i E N) is a sequence of independent processes 
such that fi as fti e i aw p s _ ; an d if X{t) denotes the decreasing rearrangement of the family 
(x!ti\t) :i,n G Nj, then X = (X(t),t > 0) is a version of X started from the configuration s. 

The self-similarity and branching assumptions (i-ii) enable us to focus on the special case 
when the fragmentation starts from a single fragment with unit size, without losing generality. 
Thus, in the sequel, we whall implicitly work under the probability measure P := Pi. 
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It is intuitively obvious that the behavior of a self-similar fragmentation should depend 
crucially on the sign of the index of self-similarity. Roughly, since fragments get smaller as 
time passes, the rate of dislocations decreases when the index is positive, whereas it increases 
when the index is negative. In the critical case when a — 0, the fragmentation process is called 
homogeneous, as the rates at which fragments split do not depend of their sizes. Homogeneous 
fragmentations are of course easier to study, and more results are known for homogeneous 
fragmentations than for general ones. 

In the case when splits occur discretely, that is when each fragment remains unchanged 
for some random time and then splits, fragmentation processes can be described in terms of 
branching Markov chains in continuous time. Further, they can be endowed with a discrete 
genealogical structure, which enables us to develop their study using general techniques from 
the theory of multiplicative cascades. Plainly, this approach fails when splits occur continuously, 
that is when each fragment breaks down instantaneously. In the sixth section of this paper, we 
shall circumvent this fundamental difficulty by performing a spatial discretization. Specifically, 
Kingman's theory of exchangeable random partitions provides the right framework for studying 
a fundamental special class of fragmentation processes, called exchangeable. This approach 
enables us to reveal the fine structure of exchangeable fragmentation and to extend results 
which were proven in the discrete setting to this general setting. Finally, the last section of this 
paper is devoted to a rather informal discussion of the duality relation involving time-reversal, 
which exists between certain pairs of fragmentation and coalescent processes. 

2 Fragmentations as branching Markov chains 
2.1 Construction of fragmentation chains 

A self-similar fragmentation process is called a (self-similar fragmentation) chain if the first 
dislocation time 

T := inf {t > : X(t) ^ (1,0,...)} 

is strictly positive a.s. (recall that we implicitly assume that at the initial time there is a single 
fragment with unit size). In this situation, the Markov property forces T to have an exponential 
distribution, say with parameter c > 0. Excluding implicitly the degenerate case when c = 
(i.e. T = oo a.s.), we may and will henceforth focus on the case when c = 1 for the sake 
of simplicity, as the general case can be reduced to the former by a linear time-substitution. 
We write v for the distribution of X(T) under So v is a probability measure on S with 
z/({(l, 0, . . .)}) = 0, which we call the dislocation law of X. 

We can then think of X as the evolution of a non-interacting particle system in ]0, oo[, 
in which each fragment with positive size is viewed as a particle and different particles have 
independent evolutions. Specifically, a particle with size r > lives for an exponential time 
with parameter r a . Then it disappears and is replaced by a cloud of smaller particles, say 
(ri, . . .) G S, such that the sequence of ratios (ri/r, . . .) has the fixed law v and is independent 
of the lifetime of the particle r. This description makes sense only when r > 0; however by 
assumption, the possible children of a particle with size have all size 0. Particles with size 
play no role, and the evolution is thus well-defined in all cases. It should be intuitively 
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clear from this description that the law of the system is entirely determined by the dislocation 
measure v and the index of self-similarity a. 

Conversely, given a real number a6l and a probability measure v on S which fulfils mild 
conditions, it is easy to construct a self-similar fragmentation chain with index a and dislocation 
law v. Roughly, the idea is to focus first on particles with size at least e > 0, and to built the 
restricted system using the standard theory of Markov chains in continuous time. Then one 
observes a compatibility property for different values of the threshold parameter e > 0, which 
allows us to take a projective limit as e — > 0. 

More precisely, assume that i/({s e S : Si > 1}) = 0, v({l, 0, . . .)}) = 0, and 

# {i : Si > e} u(ds) < oo for every s > 0. (1) 

For every e > 0, one can then construct 1 a branching Markov chain in continuous time, 
say X( £ \ with values in the space of finite atomic measures on ]e, oo[ and governed by the 
following transitions. When the process starts from say m = Yli=i with r\ > r 2 > . . . > 
r n > e, the first jump occurs at an exponential time with parameter Y^i=i r f- The configuration 
immediately after the jump is independent of the exponential time and has the same distribution 
as the point measure obtained from m by picking an atom T{ at random with probability 
proportional to rf , and replacing it by a random family of atoms, say (r»si, . . . , r,Sfc) where 
s = (si, . . .) has the law v and k = max{j e N : TiSj > e}. It is then immediately checked from 
calculations of branching rates that for every e > i] > 0, the restriction of the point measure 
process to }e, oo[ is a version of X^. By Kolmogorov's extension theorem, this ensures 
the existence of a unique law of a process X = (X(t),t > 0) with values in the space of Radon 
point measures on ]0, oo[, such that the restriction of X to }e, oof is a version of X^ for every 
£ > 0. The natural identification of a Radon point measure with the ranked sequence of its 
atoms enables us to view X as a self-similar fragmentation chain with the desired dynamics. 

We stress that in general, despite of the terminology, X is not a continuous time Markov 
chain. Indeed, even though the lifetime of a fragment is exponentially distributed with a positive 
parameter, the process may create infinitely many fragments in a finite time (but of course, 
only finitely many of them have a size greater than e for any fixed e > 0), and the infimum of 
the lifetimes of an infinite family of fragments can be zero. 

In the special case a = 0, the lifetime of each fragment is a standard exponential vari- 
able, independently of the size of the fragment, and we say that the fragmentation chain is 
homogeneous. In this situation, there is a natural connection with branching random walks 
in continuous times (cf. Uchiyama [55] )• Specifically, consider a homogeneous fragmentation 
chain X with a dislocation law v that charges only the sub-space of sequences s = (s\, . . .) with 
Sk = for k sufficiently large. Then the process 

Z®(dx) :=£<L taJfi(t) (dz), 

where the sum is taken over the fragments with strictly positive size, is a branching random 
walk. More precisely, its branching measure is the image of v by the map x —>■ —Inx. This 

1 Condition is needed to prevent a possible explosion, as otherwise it could happen that the dynamics 
would create an infinite number of atoms of size > e in finite time. 
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elementary connection has a number of interesting consequences as it essentially reduces the 
study of the class of homogeneous fragmentations associated to a dislocation laws charging only 
finite configurations, to that of branching random walks on ]0, oo[, for which a great deal of 
results are known. 

The infinitesimal generator G of a self-similar fragmentation chain has a simple and useful 
expression for so-called additive and multiplicative functionals, which is easily obtained from 
the construction described above in terms of the jump rates. 

Proposition 1 (i) Consider a measurable function f : [0, oo[^ R with f = on some neigh- 
borhood ofO, and define an additive functional A : S — > R by 

oo 

A ( S ) = Yl S= ( S 1> S 2,---)- 

i=l 

Then for every x = (x±, . . .) G S, we have 

oo „ 

GA(x) = J2l {xi>0} x? / (A(x lS )-f(x t ))u(ds). 
i=i J s 

(ii) Consider a measurable function g : [0, oo[— >]0, oo[ with g = 1 on some neighborhood ofO, 
and define a multiplicative functional M : S — ►](), oo[ by 

oo 

i=i 

Then for every x = (x±, . . .) G S, we have 

GM(x) = Vl {lj > 0} x^ / (M(x i s)-^(x i ))Kds). 



2.2 Some analytic expressions and Malthusian hypotheses 

Throughout the rest of this section, we will be working with a self-similar fragmentation chain 
X with index a G R and dislocation law v, where v is a probability measure of S satisfying the 
conditions of the preceding section. To avoid uninteresting discussions, we shall not consider 
the situation when Sj = or 1 for all i G N, v(ds)-&.s., since the fragmentation chain can then 
be identified as a branching process in continuous time. 

For future purposes, it is convenient to introduce the notation 



p := inf |p > : J ^ s p i v{ds) < oo | 



(with the convention inf = oo), and to assume from now on that p < oo. Then we define for 
every p > p 

«(p):=jf ^l-f;^i/(ds). (2) 
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Note that our assumptions ensure that k is always a continuous strictly increasing function on 
]p, oo[; n(p—) may be finite or equal to — oo. 

We then make the fundamental : 

Malthusian Hypotheses. There exists a (unique) solution p* > p to the equation n{p*) = 0, 
which is called the Malthusian exponent. Furthermore the integral 

/ s (fx)V> 

is finite for some p > 1 . 

We stress that the Malthusian hypotheses are quite weak. For instance, if si < 1, v(ds)-&.s. 
(which means that the fragments resulting from the dislocation of a particle are always strictly 
smaller than the initial particle), then we have by dominated convergence that lim^oo n{p) = 1. 
If moreover we can find some p > p such that < (this occurs for instance whenever 

f s #(s)v(ds) g]1,oo[ where #(s) := Card{i : > 0}, then the fact that k is s a continuous 
and strictly increasing ensures the existence of the Malthusian parameter. 

Throughout the rest of this work, the Malthusian hypotheses will always be taken for granted. 
The function k and the Malthusian exponent p* play a crucial role in the study of the asymptotic 
behavior of fragmentation chains, as we shall see soon. 

3 Fragmentation chains and multiplicative cascades 

In this section, we point at a representation of a fragmentation chain as an infinite tree with 
random marks. This representation can be viewed as a different parametrization of the process, 
in which the natural time is replaced by the generation of the different particles. 

3.1 Genealogical coding of fragmentation chains 

We start by introducing some notation. We consider the infinite tree 

oo 

U := (J N n , 

n=0 

with the convention N° = {0}. In the sequel U will often be referred to as the genealogical tree] 
its elements are called nodes (or sometimes also individuals) and the distinguished node the 
root. Nodes will be used to label the particles produced by a fragmentation chain. For each 
u = (ui, . . . ,u n ) G U, we call generation of u and write \u\ = n, with the obvious convention 
|0| = 0. When n > 1 and u = (u±, . . . ,u n ), we call u— := (ui, . . . , w n -i) the parent of u. 
Similarly, for every iGNwe write ui — (ui, . . . ,u n ,i) G N n+1 for the i-th child of u. Finally, 
we call mark any map from U to some (measurable) set. 

Now, we associate to each trajectory of the fragmentation chain a mark on the infinite tree 
U. The mark at a node u is the triple (£ u , a u , ( u ) where £ u is the size, a u the birth-time and 
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( u the lifetime of the particle with label u. More precisely, the initial particle corresponds to 
the ancestor of the tree U, and the mark at is the triple (1, 0, Co) where (,% is the lifetime 
of the initial particle; in particular has the standard exponential law. The nodes at the first 
generation are used as the labels of the particles arising from the first split. Again, the mark 
associated to each node % G N 1 at the first generation, is the triple cij, Q), where £j is the size 
of the i-th child of the ancestor, = + (the birth-time of a child particle coincides with 
the death-time of the parent), and Q stands for the lifetime of the i-th child. And we iterate 
the same construction with each particle at each generation. 

Clearly, the description of the dynamics of fragmentation entails that its genealogical coding 
also enjoys the branching property. Specifically, we have the following recursive description : 

Proposition 2 There exists two independent families of i.i.d. variables indexed by the nodes 

of the genealogical tree, (^ u ,,u G U^j and (e U) .,u G U), where each := (i Ul ,...,u n ,i,i G N) is 

distributed according to the law v, and each e U) . = (e ui ,i G N) is a sequence of i.i.d. standard 
exponential variables, and such that the following holds: 

Given the marks a v , ( v ), \v\ < n) of the first n generations, the marks at generation n + 1 
can be expressed in the form 



where u = (ui, . . . , u n ) and ui = (ui, . . . , u n , i) is the i-th child of u. 

Proposition |21 shows that the sizes at nodes u G U) define a so-called multiplicative 
cascade; see the pioneer works of Mandelbrot [13], Kahane and Peyriere [36J, Mauldin and 
Williams jH]; see also Liu jH] for further references. Although this multiplicative cascade 
alone does not enable us to recover the fragmentation chain (we also need the information on 
birth-times and lifetimes as it is amplified in Proposition |3J below), classical notions and results 
in this field have a key role in the study of fragmentation chains. It should be intuitively obvious 
that one can express the fragmentation chain at time t in terms of the particles which are alive 
at time t, i.e. which are born at or before t and die after t. 

Proposition 3 With probability one, for every t > and every measurable function f : 
[0, oo[— > [0, oof with /(0) = 0, there is the identity 



Proof: We have to check that all the fragments with positive size which are present at time t 
have a finite generation, i.e. result from finitely many dislocations of the initial particle. In this 
direction, let us fix some arbitrarily small e > 0, and consider the threshold operator (p £ which 
consists of removing all the fragments with size less than or equal to e. Recall from Section 2.1 
that (p £ (X) is a Markov chain, in particular the number of jumps accomplished by this chain 
before time t is finite a.s. This number obviously is an upperbound for the generation of all 
fragments with size greater than e. □ 




), 



oo 
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3.2 Intrinsic martingale 

The purpose of this section is to introduce the so-called intrinsic martingale which is naturally 
induced by Malthusian hypotheses and the genealogical coding of fragmentations, and plays a 
crucial role in the asymptotic behavior of the latter. 

Theorem 1 The process 

M n := £ , ne Z+ 

\u\=n 

is a martingale which is bounded in L 1 (P), and in particular, uniformly integrable. Moreover, 
the terminal value .Moo is strictly positive a.s. whenever v(s\ — 0) = (i.e. a fragment may 
never disappear entirely after a dislocation) . 

In the sequel, (A4 n , n G Z + ) will be referred to as the intrinsic martingale. Observe that in 
the important case when dislocations are conservative, in the sense that Yli^i s i = 1> ^(tfe)-a.s., 
then p* — 1 and A4 n = 1 for all n G Z + , and the statement is trivial. 

In general the distribution of the terminal value of the intrinsic martingale is not known 
explicitly. However, it is straightforward from the branching property that there is the identity 
in law 

oo 

Mx^Eef-MW (3) 

where £ = G N) has the law u, and A4^ are independent copies of .Moo, also independent 
of £. It is known that under fairly general conditions, such equation characterizes the law .Moo 
uniquely, see e.g. [J21E2]- We also refer to Liu [H] for information of the tail behavior of the 
solution. 

4 Some applications 

In this section, we present some results on the behavior of self-similar fragmentation chains 
which can be derived by the combination of standard techniques from the theories of Markov 
chain in continuous time and multiplicative cascades. 

4.1 Some strong limit theorems (a > 0) 

The Malthusian parameter p* and the terminal value .Moo of the intrinsic martingale have a 
crucial role in the study of the asymptotic behavior of additive functionals of the fragmentation, 
i.e. of the type 

oo 

A(x(t)) = J2xf(t)f(Mt),t), t>0 

1=1 

for some measurable function / : M? + — ► 1R + . Indeed, let us denote by Ai(t) this quantity for 
/ = 1, so Ai(t) can be viewed as the analog of the intrinsic martingale when the parameter 
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is time instead of generation. It is easy to check that for a > 0, A4(t) converges a.s. and in 
L 1 (P) to the terminal value Ai^oo of the intrinsic martingale. More generally, the branching 
property and a variation of the law of large number enable to show that under appropriate 
hypotheses on the function /, it holds that 

lim A(X(t)) = ciftMoo , in L l (V). (4) 

t— >oo 

Here, c(f) is deterministic factor which can be determined by first moment calculations. The 
terminal value of the intrinsic martingale .Moo thus appears as a kind of universal random 
weight. We refer e.g. to [H3JEIJ] and references therein for many results in this vein for certain 
branching processes. In order to avoid some technical discussion related to periodicity, we shall 
often make the assumption that the dislocation law v is non- geometric , in the sense that there 
exists no real number r g]0, 1] such that Sj G {r n ,n G Z + } U {0} for all % G N, u(ds)-a.s. 

Let us first consider the case of homogeneous fragmentation chains, i.e. with index of self- 
similarity a = 0. In this situation, each particle has a standard exponential lifetime, so infor- 
mally when n G N is large, particles at generation n are alive at times close to n (by the law 
of large numbers), and for the same reason, when t is large, the generation of particles alive 
at time t is close to [t]. In this direction, one naturally expects that the strong limit theorems 
for multiplicative cascades can be shifted without significant modifications to homogeneous 
fragmentation chains. This is indeed the shown in the following result which can be 

traced back (in a simpler setting) to Kolmogorov [3H] ; see also [2J |2H ESj for closely related 
statements. 

Proposition 4 Let f : R — > K be a continuous bounded function. In the homogenenous case 
a = 0, the following limits hold in L l (¥): 

OO 

lim^Xf (t) fir 1 In X t (t)) = M^fi-K'if)) 

t— >oo z — ' 

i=l 

and 

OO 

limY^XfVfir^ilnXM+K'ipy)) = M oo E(f(M(0,-K"(p*))), 

t— >oo z — ' 

i=l 

where jV(0, —k"(p*)) denotes a centered Gaussian variable with variance —k"(j>*). 

Roughly, the first part of Proposition 0] claims that in a homogenenous fragmentation, 'most' 
fragments decay exponentially fast with rate n'{p*). The second part is sharper and specifies 
the pathwise fluctuations. 

When the index of self-similarity a is positive, large fragments split faster than small ones, 
so the rate of fragmentation decays as time passes and one may expect a homogenisation 
phenomena. Theorem |21 below was first established by Filippov |2H] in the special case when 
the dislocation measure v is conservative, i.e. when *Y^i=\ Si = 1 f° r ^" a - e - s ( see a l so Brennan 
and Durrett Recall that in this case, p* = 1 and A4 n = 1. The general case of non- 

conservative dislocation measures was recently proved by Bertoin and Gnedin j!5j . 
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Theorem 2 Suppose that a > and that the dislocation law v is non- geometric. Then for 
every bounded continuous function f : M + — ► K. 

lim Vxf (t)f(t^ a Xi(t)) = / f(y)p(dy) , zn ^(P), 

where Aioo is the terminal value of the intrinsic martingale and p is a deterministic probability 
measure. More precisely, p is determined by the moments 

y ak p(dy) = - . . ^\ — - for keN, 

]o,oo[ aiK'[p*) K{p* + «)••• «(p* + (A - l)a) 

(it/it/i the usual convention that the right-hand side above equals l/(<W(p*)) for k = 1). 

Comparison with the homogeneous case of Proposition 0] shows that the size of a typical 
fragment now decays as a power function of time. It is also interesting to observe that the limit 
is much more sensitive to the dislocation law than in the homogeneous situation: the function 
k can be recovered from the limit measure p, whereas in the homogeneous case, the exponential 
rate of decay just depends on the derivative of k at the Malthusian parameter. 

Let us just sketch an argument for the proof of Theorem |2] based on moment calculations. 
Kolmogorov's backwards equation combined with Proposition ^yields, in the special case of a 
power function, the equation 



dt 

The solution is given in the form 



( 5>f- Q (t) ) = -«(p- Q )E ( $>f(t) 

i=l / \i=l 



(oo \ oo /_,w 

8=1 / n=0 

where 7(0, p) = 1 and for n > 1 

n-1 

j(n,p) = J^K(p + ak). 



k=0 

Asymptotics can be obtained using techniques from complex analysis, and one gets for p = 
p* + ak with k G N, that 

V-^i / aK'{p*) K{p* + a) ■ ■ ■ K{p* + [k — l)a) 
This shows that 

(°° \ poo 

^Xf(t)f(&'>X i (t))j = j( f(y) P (dy) 

when f(x) = x ka , and then for any continuous function / bounded by a power function. 
Finally, /^-convergence in Theorem El follows from the meta-limit theorem (J3J). We refer to ^H] 
for details. 
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4.2 Extinction and formation of dust for a < 



Intuitively, when the index of self-similarity is negative, fragments with small sizes are subject 
to high splitting rates, and this makes them vanish entirely; see e.g. pU lE?T) l3o] f° r some works 
formalizing this intuition. 

Proposition 5 Suppose 

Card {i : Si = 1} u(ds) < 1 . 



s 

Then the following assertions hold with probability one: 

(i) For a < 0, X(t) = (0, . . .) for all sufficiently large t. 

(ii) When k(—cx) > 0, for almost every t > 

Card (jeN: Xj(t) > 0} < oo . 

We stress that in general, no matter what the value of a is, there may exist random instants 
t at which 

Card (jeN: X 3 {t) > 0} = oo . 
For instance in the case when the dislocation law fulfills 

v (xj > for all j G N) = 1 , 

then with probability one, there occur infinitely many sudden dislocations in the fragmentation 
chain X, each of which produces infinitely many terms. This does not induce any contradiction 
with Proposition El (ii) when k(— a) > 0, because informally, as the index of self-similarity is 
negative, we know that fragments with small size vanish quickly. 

It would be interesting to have information on the distribution of the extinction time 

C:= inf{t>0:X(t) = (0,0,...)}, 

however it does not seem possible to express this law in a closed form. Nonetheless, we point 
that an application of the branching property at the first dislocation yields the identity in 
distribution 

C^e + max^C,, (5) 

where e is a standard exponential variable, (£,-, j G N) is distributed according to u, G N) is 
a sequence of independent copies of (, and finally, e, G N) and (Q,j G N) are independent. 
We refer to [3j for a survey of this type of equations in distribution. 

It is interesting to study in further details the extinction phenomenon in the case when the 
dislocation measure is conservative, i.e. when 



j>^l) = 0. (6) 



j=i 
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It is easy to deduce by iteration that for every n e N, the total mass of particles at the n-th 
generation is conserved, i.e. 

^2 = 1 , a.s. 

\u\=n 

Turning our interest to the total mass of particles at time t, we introduce the quantity 

oo 

D(t) :=l-J2Mt), 

n=l 

which can be viewed as the total mass of dust, that is of infinitesimal particles at time t. One 
could be tempted to believe that the assumption (0) would entail D = 0; indeed it is easy 
to check that this holds when the index of self-similarity of the fragmentation is nonnegative. 
However Proposition El shows that for negative indices of self-similarity, D reaches 1 at a finite 
time a.s. 

Proposition 6 Suppose (OJ) holds. The following assertions hold with probability one: 

(i) D is a continuous increasing process which reaches 1 in finite time. 

(ii) If := Card{i : Xi(t) > 0} denotes the number of fragments with positive mass at time 
t, then 

POO 

/ l{#(t)<oc}dD(t) = 0. 

JO 

This statement again reflects the fact that, loosely speaking, dislocations occur faster and 
faster as time passes. Observe that it entails that almost-surely, there exists uncountably many 
times at which there are infinitely many fragments with positive size, which may be rather 
surprising (for instance in the case when dislocations are binary, i.e. produce exactly two 
fragments; see also Proposition 

Recall that when the dislocation law v is conservative, n(p) > for every p > 1, and 
Proposition El shows that for a < —1, at each fixed time t there is only a finite number of 
fragments with positive size in the system. Combining this observation with Proposition El 
entails that the random measure dD(t) is singular with respect to Lebesgue measure on R + . 
Haas [32| established the following sharper result. 

Proposition 7 Suppose (OJ) holds and that there exists k G N such that z/(s/c +1 > 0) = (i.e. 
each dislocation splits a fragment into at most k pieces). The following assertions hold with 
probability one: 

(i) If a e] — 1,0[ ; then D is absolutely continuous. 

(ii) If a < —1, then the random measure dD{t) is supported by a set with Hausdorff dimension 
l/\a\. 

5 Fragmentation chains and general branching processes 

We now turn our attention to a slightly different aspect of fragmentation chains, by pointing 
at a connection with general (i.e. Crump-Mode- Jagers) branching processes; see |34[ HU] and 
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references therein. Specifically, we use the genealogical coding of the fragmentation chain of 
Section 2, and for each node u EU, we think of a u := — ln£ u as the birth-time of an individual. 
Every individual which is born has an infinite lifetime, and gives birth to children according 
to a random point process More precisely, if si > S2 > . . . denote the sequence ranked 
in the decreasing order of the sizes of the fragments that result from the first dislocation of 
a fragment with size s (recall that we assume that s > s±, i.e. dislocations always produce 
fragments of smaller sizes), the nth child of that individual is born at time — lns n . In other 
words, each fragment with size s is viewed as an individual u born at time a u = — In s and 
which has offspring described by the point process 

Vu ([0, t}) = Card {n G N : s n > e - *} . 

The particle system obtained in this way is then a general branching process with reproduc- 
tion intensity 

„ 00 

//(*) = E( m ([Q,t])) = / Y)l K >e-^(<fe), t>0. 

J s n=1 

This point of view is useful to investigate problems related to situations when fragments with 
size less than a certain fixed parameter e > are instantaneously frozen. The fragmentation 
process then terminates when the system only constists in particles with size less than e. This 
setting arises for instance in the mining industry where fragmentation is needed to reduce rocks 
into sufficiently small particles. For this purpose, rocks are broken in crushers and mills by 
a repetitive mechanism. Particles are screened so that when they become smaller than the 
diameter of the mesh of a thin grid, they are removed from the process. 

Imagine the instantaneous dislocation of a block of size s into a set of smaller blocks of 
sizes (si, S2, ...) requires an energy of the form s^ip^si/s, S2/S, . . .), where if : S — > R is a cost 
function and (3 > a fixed parameter. We are interested in the total energy-cost of the process 
that finishes when all particles have size less than e > 0: 

S{e) := X>fe>^M«.), 

where £ u . = {^u,n/iv,i n £ N) is the sequence of ratios of the sizes of the fragments resulting 
from the dislocation of the fragment labelled by the node u and the size of that fragment (see 
Proposition |2J) . 

The following asymptotic result for the energy cost as e — > has been derived recently by 
Bertoin and Martinez [TS] from the work of Nerman ^1 o n strong limit theorems for general 
branching processes. 



Theorem 3 Suppose that the dislocation law is non-geometric and fulfils J s \(p(s)\u(ds) < 00. 
Then lim £ _^ £(e) := £(0+) in L X (P) for (3 > p* , whilst for f3 < p* it holds that 

\ime pt - l3 S(e) = M J° . / <p(s)v(d&) m L\F) . 

e-0 K ' (p* - /3)K'(p*) J s ^ K ' V ' K ' 
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We may also consider the distribution of the small particles that can go across the grid. 
Specifically, we would like to get information about the random finite measure on ]0, 1[ 

g £ (dx) := ^ ^.-^J^^/eW) 

which can be viewed as a weighted version of the empirical measure of the particles taken at the 
instant when they become smaller than e and then renormalized. In the setting of the general 
branching process associated to the fragmentation chain, we can re-express g £ in the form 

oo 

(&, /> = ^^e-^«»/(e i --)l {ffu< K ff „„ } , 

u£U n=l 

where t — — Ins and / :]0, 1[— > M + denotes a generic measurable function. 

Proposition 8 Suppose that the dislocation law v is non- geometric. As e — > 0, g £ converges 
in probability to Ai^n, where g is a deterministic probability measure on [0, 1] given by 




6 Fragmentations with instantaneous dislocations 

In this section, we turn our attention to the situation when dislocations can occur instan- 
taneously. Examples arise naturally e.g. in the study of Brownian motion and Continuum 
Random Trees; see jU El E31 E3 E|- Then one cannot consider the first dislocation of a 
particle, and the genealogical structure of the process is no longer discrete, which impedes 
the representation as a multiplicative cascade. Our purpose is two-fold; first we would like to 
characterize such fragmentation processes and their structure, and second we aim at extending 
properties which have been established for fragmentation chains to this more general setting. 

We shall focus on the case when the process X takes values in the space of mass-partitions 

A mass-partition s £ V m can be thought of as the ranked sequence of the masses of the fragments 
of some object with unit total mass; the case when s i < 1 corresponds to the existence 
of dust, i.e. a part of the object with positive mass has been reduced to infinitesimal particles 
each with zero mass. 

It might be natural to try to investigate fragmentations with instantaneous dislocations as 
limit of fragmentation chains when the parameter of the exponential lifetime of a particle with 
unit mass tends to oo (so the lifetime of that particle tends to 0). However this approach is 
far from being easy, and we shall rather follow a different route initiated by Kingman [37 and 
then further developed by Pitman [HU] in the framework of coalescents (see also Schweinsberg 
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6.1 Exchangeable partitions and interval representations 

Kingman's idea is to encode mass-partitions by random partitions of N. Specifically, consider 
some space E endowed with a probability measure g, and disjoint measurable subsets Ei, . . . 
with positive ^-measure. Set E := E\ (U k >iE k ), so (E ,Ei, . . .) is a partition of E. The sets 
Ei, . . . are viewed as fragments of E and E as the set of dust. 

Then let Ui, . . . be a sequence of i.i.d. variables in E with law g, and define a random 
partition 7r of N as follows. If Ui G E , then {i} is a singleton of the partition 7T, and the 
other blocks of n are of the type n k = {£ E N : Ue <E E k } for A; > 1. The strong law of large 
numbers enables us to recover the sequence of the masses of the partition of E as the asymptotic 
frequences of the blocks of ir; more precisely 

g(E k ) = lim -Card {£ e ir k : £ < n} . 

We stress that only a special class of partitions of N can arise in this framework, as blocks of 
such partitions always have asymptotic frequences. More importantly, the random partitions 
resulting from Kingman's construction are exchangeable, that is their distribution is invariant 
by the action of permutations. Indeed the action of a permutation a amounts to a permutation 
of the indices in the sequence of i.i.d. variables, and since U a n), U a ( 2 ), ■ ■ ■ has the same law as 
Ui, U2, ■ ■ ., we see that the image of tc by the action of a has the same law as n. 

Kingman [37| (see also Aldous pQ for a simpler proof) has shown that for any random 
exchangeable partition 7 of N, the blocks of 7 possess asymptotic frequences, and more precisely, 
7 has the same distribution as some partition 7r constructed as above for a certain random 
partition of a space E. Specifically, one can take E =]0, 1[ endowed with the Lebesgue measure, 
let Ei, ... be the interval components of some random open set G C]0, 1[ and E =]0, 1[\G. 
Then Ui, . . . is a sequence of i.i.d. uniform variables which is independent of the random open 
set G. 

Now consider a fragmentation process X in the sense of Definition Q and assume that X 
takes values in the space of mass-partitions. For our purpose, it is convenient to think of X 
in terms of a fragmentation of the unit interval, in the sense that there is a Markov process 
(G(t),t > 0) with values in the space of open sets in ]0, 1[, such that G(t) C G(s) when s < t 
and for each t > 0, X(t) is the ranked sequence of the interval components of G(t). The fact 
that such a representation exists is explained in [7| I12j. 

Next, let Ui, ... be a sequence of i.i.d. uniform variables which is independent of (G(t),t > 0), 
and for each t > 0, write Tl(t) for the random exchangeable partition of N such that two distinct 
indices i,j belong to the same block of U.(t) if and only if Ui and Uj belong to the same interval 
component of G(t). We make the key observation that for every t > 0, given an interval 
component, say I of G(t), if B = {i G N : Ui G 1} denotes the block of U(t) corresponding to /, 
then conditionally on B, (U,i G B) is a sequence of i.i.d. uniform variables in /. Essentially, 
this observation implies that the partition- valued process II = (U(t),t > 0) is Markovian. Thus 
to each fragmentation process X we can associate a Markov process IT with values in the space 
of partitions of N such that IT(t) gets finer as t increases. More precisely, the law of II is 
invariant by the action of permutations and the branching property of X is transferred to IT. 
Thus II can be thought of as a fragmentation process with values in the space of partitions of 
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N. 

Conversally, given an arbitrary fragmentation process II with values in the space of partitions 
of N and which is invariant by the action of permutations, Kingman's theorem enables us to 
recover a fragmentation process X with values in V m by considering the asymptotic frequences 
of the blocs of II(t). We refer to [7j for a precise argument. 

6.2 The structure of homogeneous fragmentation processes 

In this section, we suppose that the index of self-similarity is a = 0. In this situation, the 
fragmentation process II with values in the space of partitions of N associated to X as above 
enjoys a crucial additional property : for every n G N, the restricted process II|[ n ] is Markovian, 
where for a partition tt of N, 7rir n ] denotes the restriction of 7r to [n] := {1, . . . ,n}. Informally, 
this follows from exchangeability and the fact that the dislocation rates of blocks of IT do 
not depend on their sizes. Since the space of partitions of [n] is finite, II|[ n ] is a continuous 
time Markov chain, whose evolution is thus specified by its jump rates. Combining this with 
Kingman's theory of exchangeable random partitions enables us to reveal the structure of ho- 
mogeneous fragmentation processes; see [TT] for details. Roughly, one gets that homogeneous 
fragmentations result from the combination of two different phenomena: a continuous erosion 
and sudden dislocations. The erosion is a deterministic mechanism, analogous to the drift for 
subordinators, whereas the dislocations occur randomly according to some Poisson random 
measure, and can be viewed as the jump-component of the fragmentation. 

More precisely, let us first focus on dislocations. We call dislocation measure a measure v on 
V m , which gives no mass to (1,0,...) and fulfills 

/ (l-x 1 )v(dx) < oo. (7) 

Then, consider a Poisson random measure on V m x N x M + , 

oo 

1=1 

with intensity v <S> # <S> dt, where # denotes the counting measure on N. One can construct 
a pure jump process (Y(t),t > 0) in V m which jumps only at times t(i) at which there is an 
atom (A(i),k(i),t(i)) of the Poisson measure. More precisely, the jump (i.e. the dislocation) 
induced by such an atom can be described as follows. 

The mass-partition Y(t(i)) at time t(i) is obtained from that immediately before t(i), i.e. 
Y(t(i)—), by replacing its k(i)-th term, viz. Y^(t(i)—), by the sequence Yk^(t(i)—)A(i), and 
ranking all the terms in the decreasing order. For instance, if 

Y(t(i)-) = Q, i ^ 0, . . , fc(i) = 2 and A(i) = 0, . . 

then we look at the 2-nd largest term in the sequence Y(t(i)— ), which is j, and split it according 
to A(i). This produces two fragments of size ^ and ^g, and thus 

Y(t(i)) = (-—,—,— ,0,... 1 . 
v v ;; V 3 16 12 16 / 
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Next, call erosion coefficient an arbitrary real number c > 0, and set X(t) = e~ ct Y(t). So 
X is obtained from Y by letting the fragments of the latter be eroded at constant rate c. 
The process (X(t),t > 0) is again a homogeneous fragmentation, and conversely any homoge- 
neous fragmentation process can be constructed like this. In conclusion, the distribution of a 
homogeneous fragmentation is entirely specified by its erosion rate and its dislocation measure. 

In the special case when c = and v is a finite measure (and, in particular, a probability), it 
is easy to check that X is then a homogeneous fragmentation chain with dislocation law v in the 
sense of Section 2. Note however that the condition (JJJ) allows v to be infinite, which corresponds 
to the situation when particles dislocate instantaneously. Informally, mass-partitions s for 
which 1 — s\ is small should be thought of as small, in the sense that a small mass-partition 
produces one large fragment and all the remaining ones are small (in particular, the mass- 
partition (1, 0, . . .) has to be viewed as a neutral element for dislocations). So, roughly speaking, 
condition (JJJ) allows infinitely many small dislocations, but guaranties that the accumulation of 
these small dislocations does not reduce instantaneously the initial mass into dust. This bears 
obvious similarities with subordinators, which are constructed by the Ito-Levy decomposition 
from atoms of a certain Poisson random measure on R + , see Chapter 1 in Pj. The intensity of 
this Poisson random measure is given by the so-called Levy measure A of the subordinator, and 
the integrability condition J R (1 A x)A(dx) < oo for a measure on M + to be the Levy measure 
of some subordinator is the necessary and sufficient condition for the summability of the atoms. 

The Poissonian structure of homogeneous fragmentation is a fundamental tool which enables 
to circumvent difficulties related to the absence of a discrete genealogical structure. However, 
although the law of a homogeneous fragmentation process is characterized by its erosion rate and 
dislocation measure, in general we do not know how to describe explicitly e.g. the distribution 
of the process at a fixed time. The next crucial tool for the study lies in the fact that partial 
but most useful information can be derived from a so-called size-biased sampling, and it turns 
out that the law of the latter is simple to formulate. 

Recall that we may represent X in terms of fragmentation of the unit interval. So we consider 
a Markovian family (G(t),t > 0) of nested open subsets of the unit interval, in particular for 
every t > 0, X(t) is the ranked sequence of the lengths of the interval components of G(t). Now 
suppose that U is a uniform random variable on ]0, 1[, which is independent of (G(t),t > 0), 
and for every t > 0, denote by x{t) the length of the interval component of G(t) which contains 
U (with the convention that x(f) — if U £ G(t)). In other words, the process x — (x(0> ^ — 0) 
gives the size of the fragment containing a point which has been tagged independently of the 
fragmentation process; one refers to x as the process of the tagged fragment. Note that x(f) is 
a size-biased pick from the sequence X(t) = (X 1 (t), . . .), i.e. there is the identity in law 

X {t) = X K (t) , 

where K is an integer valued variable whose conditional distribution given X(t) is 

¥(K = k | X(t)) = X k (f) , k = l,.... 

It turns out that the process of the tagged fragment is closely related to a subordinator (i.e. 
an increasing process with independent and stationary increments; see jH] for background) that 
can be characterized explicitly in terms of the dislocation measure and the erosion coefficient. 
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Theorem 4 The process 

a{t) = -lnx(t), t>0 

is a subordinator, whose law is determined by its Laplace transform of its one- dimensional 
distributions. We have 

E ( X (t) q ) = E(exp(-ga(t))) = exp (-t«(g + 1)) , t, q > , 

where the function k is given in terms of the erosion rate c and the dislocation measure v by 
the identity 



K(q) := cq + j^ (l - «A i/(ds) 



g > l 



Many features (like asymptotic behavior) of fragmentation process can be read on properties 
of tagged fragments, the combination of Theorem 0] and the theory of subordinators provides 
the key to many results on homogeneous fragmentations. 



6.3 Additive martingales and applications 

The independence and stationarity of the increments of subordinators entail that for every 
q > 0, the process 

exp(-gor(t) + K(q + l)t) = exp(t/t(g + l))x q (t) , t>0 

is a martingale. As the tagged fragment x(t) is picked at random from the mass-partition X(t) 
by size-biased sampling, it follows (take q = p + 1) that 

M(p,t) := exp(tK(p))J2x?(-t) > * > 

i=l 

is a nonnegative martingale for every p > p (this can also be checked directly from the branching 
and scaling properties). 

One refers to M(p,t) as an additive martingale. Plainly M(p,t) converges a.s. as t — ► oo, 
and in order to investigate the asymptotic behavior of homogeneous fragmentation processes, 
it is important to know if the limit of this martingale is strictly positive or zero. We shall 
investigate this question in the special case when the dislocation measure is conservative and 
there is no erosion, i.e. we assume from now on that 

oo 

c = and s« = 1 v(ds)-&.e. (9) 

i=i 

Observe that the Malthusian hypotheses are then automatically fulfilled and the Malthusian 
parameter is p* = 1. A first step in the analysis is the following elementary lemma. 
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Lemma 1 Assume (GJ). The function p — > k(p)/p reaches its maximum at a unique location 
p > 1, which is the unique solution to the equation 

pn'(p) = K,(p) . 

More precisely, the function p — > n(p)/p increases on ]p,p[ and decreases on ]p,oo[, and the 
value of its maximum is n'(p) = n{p)/p- 



Proof: We first point out that the function k is concave and increasing. It follows that 

the function p — > pn'(p) — n(p) decreases on ]p, oof. (10) 

Indeed, this function has derivative pk"(p), which is negative since k is concave. Recall from 
(JHJ) and (jHJ) that k(1) = 0; on the other hand, it is obvious that lim^oo K(q)/q = 0, hence the 
function p — > n{p)/p has the same limit at 1 and at oo, so it reaches its overall maximum at 
a unique point p > 1. In particular, we deduce from (fTUj) that the derivative of p — > n(p)/p is 
positive on }p,p[ and negative on ]p, oof. Finally, the derivative must be zero at p, which entails 
that the overall maximum is given by n'(p) = n(p)/p. □ 

We may now state the main result of this section which can be proved using the Poissonian 
structure of homogeneous fragmentations and stochastic calculus, see [12] for details. 

Theorem 5 Assume (0J). For every p (=]p,p[, the martingale M(p, •) is bounded in L X (P) and 
its terminal value is strictly positive. 



Just as Theorem^ for the intrinsic martingale, Theorem has crucial role in the study of 
the asymptotic behavior of homogeneous fragmentation. First, we specify the rate of decay of 
the largest fragment, refering to Biggins for a similar result in the framework of branching 
random walks. 



Corollary 1 Assume (GJ). It holds with probability one that 

lim-lnX^t) = -«'(p) = 

t^oo t p 

Proof: For every p > p, we have 

oo 

exp(tK(p))X>(t) < exv(tK(p))J2xm 

and the right-hand side remains bounded as t tends to infinity. Hence 

limsup-lnXxO) < 

t^oo t p 



and optimizing over p yields 



limsup-lnXi(t) < -^M 

t^oo t p 
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On the other hand, for every p &]p,p[ and e > sufficiently small, we have the lower bound 

oo oo 

exp(t«(p)) < X?(t) exp(t«(p)) ^Xf- £ (*) • 

i=l i=l 

We know that both limits 

oo oo 

lim exp(tn(p)) } Xf(t) and lim exp(t/t(p — e)) N Xf _e (t) 

t— >oo ^— ' t^oo ' ^ 

i=l »=1 

are finite and strictly positive a.s., and we deduce that 

J, , i — k(» — e) 

liminf -InX^t) > 



t—>oo t £ 

We take the limit of the right-hand side as e — > 0+ and then as p tends to p to conclude that 

lim inf - In X\ (t) > —k'(j>). 

t—*oo t 

Now, this quantity coincides with —n(p)/p, as we know from Lemma ^ □ 

It is interesting to compare Corollary ^ with Proposition E[i), which claims that the size of 
most fragments decays exponentially fast with rate k'(p*) = The size of largest fragment 

thus also decays exponentially fast, but with a slower rate k'(j>) < We refer to Berestycki 

8j for the multi-fractal analysis of the exponential rates of decay of fragments in homogeneous 
fragmentations. 

Alternatively, one can also establish Theorem El by discretization, using the following con- 
nection with branching random walks and classical results on the latter; see ^HI- If we consider 
the point measure with atoms at the logarithms of the fragments 



1=1 



we can think of the discrete skeleton {Z^ n \n G N) as a non-interacting particle system. Specif- 
ically, particles evolve independently one of the other, and at each step, each particle, say y, is 
replaced birth to a random cloud of particles y+Z, where the law of Z is that of (In Xi(l),i G N). 
This means that (Z( n \n > 0) is a branching random walk in the sense of [5] , [22] > ■ • ■ The latter 
have been throughoutly studied in the literature, and many of their properties can be translated 
to homogeneous fragmentation processes. For instance, we can derive precise information on 
almost sure large deviations for the empirical measure; using a genuine result of Biggins [21] 
for branching random walks. 

Corollary 2 Assume (0) and that the dislocation measure v is non- geometric, and for p G 
}p,p[, let M(p, oo) denote the terminal value of the uniformly integrable martingale M(p, •). // 
/ : R — > R is a function with compact support which is directly Riemann integrable, then 

lim Vie-^-^ f f(tK>(p)+y)ZU{dy) = f°° /(y)e™dy. 

uniformly for p in compact subsets of ]p, p[, almost surely. 
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In a different direction, one can also use the discretization techniques to estimate the prob- 
ability of presence of abnormally large fragments as time goes to infinity; see JH] for details. 

Corollary 3 Assume (GJ) and that the dislocation measure v is non- geometric, fix two real 
numbers a < b and take p > p. Then as t — ► oo 

P (3i E N : Xi(t) E [ae- tK ' {p) ,be- tK '^ ~ Kt -i^ e t{ P n'{j>)-K{p)) ^ 

where K is some positive and finite constant depending on a, b and the characteristics of the 
fragmentation. 



6.4 Changing the index of self-similarity 

So far, we have only been able to study fragmentations with instantaneous dislocations in the 
homogeneous case, i.e. when the index of self-similarity is a = 0. In this section, we present a 
simple transformation that changes a homogeneous fragmentation X into a self-similar one with 
an arbritrary index of self-similarity, X^ a \ completing the construction of general self-similar 
fragmentation processes. We refer to ^2] for details. 

In this direction, it is convenient to start from an interval-representation (G(t),t > 0) of 
some homogeneous fragmentation X as in section 6.1. For every y e]0, 1[, let I y (t) denote the 
interval component of G(t) that contains y if y G G(t), and I y (t) = otherwise. We write |/| 
for the length of an interval / C]0, 1[, and for every y g]0, 1[ we consider the time-substitution 

T { y a \t) : = inf ju > : jf \I v (v)\- a dv > tj . 

Because the open sets G(t) are nested, we see that for every y, z e]0, 1[, the intervals I y {T^ a) (t)) 

and I z (T^ a \t)) are either identical or disjoint, so the family | I y (Ty a \t)), < y < l| can be 

viewed as the interval components of an open set G^ a \t). It is straightforward that the family 
(G^ a \t),t > 0) is nested. More precisely, if we write X^ a \t) for the ordered sequence of the 
lengths of the interval components of G'"'(i), then (X^ a \t), t > 0) is a self-similar fragmentation 
with index a. 

Any self-similar fragmentation X^ can be constructed from some homogeneous one X as 
above, and this construction can be inverted. In particular, the distribution of X^ is entirely 
determined by the index of self-similarity a, and the erosion coefficient c > and the dislocation 
measure v of the homogeneous fragmentation X. 

The key tool which is needed to extend the results of Sections 3-5 to self-similar fragmentation 
with instantaneous dislocations is provided by the stochastic structure of the process of the 
tagged fragment, {t) ,t > 0). Recall that in the homogeneous case, (x{t),t > 0) can 
be described as the exponential of a subordinator. The construction above of a self-similar 
fragmentation from a homogeneous one by time substitution enables us to derive the law of the 
tagged fragment X {t) f° r a self-similar fragmentation. 
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Specifically, let a = (a(t),t > 0) be a subordinate* with Laplace exponent k(1 + •) given by 
(jHJ). Introduce the time-change 

r(t) = inf |m : J exp(acr(r))dr > t j- , t > 0, 

and set ( t = ex v(~°'( T (J : ))) (with the convention that Q = if r(t) = oo). Then the processes 
(Cut > 0) and (x {a) {t),t > 0) have the same law. In particular, this shows that the process 
of the tagged fragment is a decreasing self-similar Markov process as introduced by Lamperti 
|40j . see the survey [2U| and the references therein. 

Results stated for self-similar fragmentation chains in Sections 4-5 can be extended verbatim 
to self-similar fragmentation processes; see [T31 ITB^ ITS} I32~j. 



7 Duality with certain coalescent processes 

Coagulation processes are used as models to describe the evolution of particle systems in which 
pairs (or, more generally, families) of particles merge as time passes. At this level of generality, 
a simple time-reversal provides an obvious connection with fragmentation. 

However in practice, just as for fragmentation, one has to make some restrictive assumptions 
on the dynamics of coagulation in order to deal with processes that can be studied mathemat- 
ically. First, one assumes that particles are determined by their masses (i.e. a positive real 
number such that the total mass is a preserved quantity when merging occurs). So there is no 
geometry involved in the system; physicists call such models mean-field. Second, one assumes 
that the evolution is Markovian, and third, that the rate at which a family of particles merges 
only depends on the particles in this family, and not on the other particles in the system. We 
shall call stochastic coalescent a coagulation process that fulfils these requirements, and refer 
to the survey by Aldous [2] and the references therein for much on this notion. 

The first two requirements (mean-field and Markov properties) are clearly compatible with 
time-reversal. However, even though the third one bears some vague resemblance with the 
branching property, there is in general no reason why it should yield the latter after time- 
reversal. Despite the absence of a general result on duality by time-reversal between frag- 
mentation and coalescent processes, there are nontheless several important examples for which 
duality holds. It would be very interesting to establish general critera for duality; we leave this 
question open as a challenge, and now conclude this survey by discussing some examples. 

The first example is the simplest, it can be constructed from a sequence U\, . . . of independent 
uniform variables on [0, 1] and an independent Poisson process N = (N t , t > 0). Specifically, for 
each time t > 0, write F(t) for the ranked sequence of the lenghts of the interval components 
of the random open set ]0, : 1 < i < N t }. It is easy to see that F = (F(t),t > 0) is a 

self-similar fragmentation chain with index a = 1. More precisely, its dislocation law is that 
of the random mass partition (1 — U/2, U /2, 0, . . .) where U is uniformly distributed on [0, 1]. 
Then it can be checked (see ^3] for details) that the exponential time reversal C(t) := _F(e - *) 
transforms the fragmentation process F into a coalescent process C. Specifically, when there 
are n > 2 particles in the coalescent C, the first coagulation occurs after an exponential time 
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with parameter n, and the pair of particles involved is uniformly distributed amongst the 
n{n — l)/2 possible pairs. One reckognizes a variation of the celebrated coalescent of Kingman 
[3*7] (which is used to describe the genealogy of large populations and has many applications 
in Biology). More generally, a similar duality holds for self-similar fragmentation chains with 
index a = 1 and dislocation law given in terms of certain n-dimensional Dirichlet distributions. 
These fragmentation and coalescent processes appear in the genealogy of Yule processes; see 

Eg. 

The second example concerns the additive coalescent, a coagulation process which arises 
e.g. as a model for the formation of drops of rain in clouds. Roughly speaking, in an additive 
coalescent, any pair of particles, say (x, y), merges at rate x+y, independently of the other pairs 
in the system. Evans and Pitman [2B] have observed that if C +,n = (C +,n (t),t > 0) denotes the 
process started from the monodisperse initial condition which consists in n particles each with 
mass 1/n, then as n — > oo, the process (C +,n (t + |lnn),t > — |lnn) converges in the sense of 
finite dimensional distributions to C + := (C + (t),t G R). The latter is known as the standard 
additive coalescent. Again the exponential time-change, F(t) = C + (— hit), transforms the 
coalescent into a fragmentation process which is self-similar with index a = 1/2. We refer to 
jH QUI for more on this topic. In this direction, we also mention that Miermont j3H] has shown 
that the exponential time-change also transforms some (but not all) non-standard additive 
coalescent into fragmentation processes; however the latter are not self-similar in general. 

The final example is due to Pitman j3U] who established a remarkable duality between certain 
coagulation and fragmentation operators based on adequat Poisson-Dirichlet distributions. An 
important special case involves Ruelle cascades and the Bolthausen-Sznitman coalescent. The 
former have been introduced in as a tool for studying Derrida's Generalized Random 
Energy Model of spin glass; and the latter in to describe the time-reversed dynamics in 
Ruelle's cascades. The Bolthausen-Sznitman coalescent also appears naturally in the genealogy 
of Neveu's branching process, see [T7]. Pitman's duality result enables us to view Ruelle's 
cascades as a time-inhomogeneous fragmentation process, see also Basdevant jH] for a recent 
development in this direction. 
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